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We determine, by means of density functional theory, the stability and the structure of graphene- 
nanoribbon (GNR) edges in presence of molecules such as oxygen, water, ammonia, and carbon 
dioxide. As in the case of hydrogen-terminated nanoribbons, we find that the most stable arm- 
chair and zigzag configurations are characterized by a non-metallic/non-magnetic nature, and are 
compatible with Clar's sextet rules, well known in organic chemistry. In particular, we predict 
that, at thermodynamic equilibrium, neutral GNRs in oxygen-rich atmosphere should preferentially 
be along the armchair direction, while water-saturated GNRs should present zigzag edges. Our 
results promise to be particularly useful to GNRs synthesis, since the most recent and advanced 
experimental routes are most effective in water and/or ammonia-containing solutions. 

PACS numbers: 71.15.Mb, 71.20.Tx, 73.20.At 



I. INTRODUCTION 

Narrow graphene nanoribbons^' are carbon allotropes 
that promise to combine the high carrier mobility of 
graphene^ with a semiconducting nature due to quan- 
tum confinement in the lateral direction. The resulting 
systems are likely to be suitable for direct application 
as channel in field-effect transistors and other semicon- 
ducting devices. The design of GNRs, and more gener- 
ally of graphene samples with controlled edges, is cru- 
cial in determining their electronic and transport prop- 
erties. Among the variety of reported experimental ap- 
proaches, three routes, intensely explored in the last cou- 
ple of years, seem at the frontier of current research: i) 
lithographic patterninj^^ ii) sonochemical methods^* 
Hi) metallic nanoparticle chemical etching (unzipping) 
of graphite/graphene (nanotubes^F^. Gaining control 
of the edge profiles and structures is of course one of the 
most important issues to solve, somehow defining the po- 
tential impact of th ese approaches. The observation of 
graphene edgef^EH and/or GNEE^ES profiles at atomic- 
scale resolution indicates that the mentioned experimen- 
tal methods are capable of producing smooth, likely non- 
chiral (i.e. no mixing of zigzag and armchair direction), 
edges. However, little is known on how to achieve either 
zigzag or armchair edges, and about the actual chemical 
termination of the unsaturated carbon bonds. 

In this regard, the huge number of theoretical works 
that recently appeared in the literature on GNRs^^ZI 
has almost exclusively focused on single-hydrogen satu- 
ration of the edge carbon atoms. This was mainly due 
to the surprising discovery of spin-polarized electronic 
states localized on the edges of zigzag GNRs | 17 | 1S | , pos- 
sibly turning them into half-metals under high electric 
fields^^l, and spurring great interest in view of possi- 



ble applications in future spintronics. However, we have 
recently^SMl pointed out that the single-H-terminated 
zigzag GNRs, and thus all those fascinating electronic 
and magnetic properties they exhibit, might only be sta- 
ble at very difficult experimental conditions. 

The hydrogen case-study is presently extended to 
GNRs terminated with other relevant atoms and 
molecule^2H2SI j n fact, as mentioned above, experiments 
are often performed in aqueous or ammonia-containing 
solutions, or in oxidizing conditions. More generally wa- 
ter and oxygen contamination (for example in air) is ex- 
tremely important in determining the electronic prop- 
erties of graphitic systems. For example, the chemical 
reactions that induce the unzipping of carbon nanotubes 
are thought to leave O-terminated zigzag GNRs as end 
products^. In this work we thus report density func- 
tional theory (DFT) calculations of the energetics and 
structures of various GNR edges terminated with atmo- 
spheric molecules, such as oxygen, water, and ammonia, 
as functions of the respective chemical potentials. The 
work is organized as follows: in section [H] we describe 
the theoretical method and define the relevant thermody- 
namic quantities; in section [TTI| we present our systematic 
study on oxygen-terminated zigzag and armchair GNRs. 
In section |TV| we report our results on selected GNRs in 
presence of water and ammonia molecules; discussion on 
the electronic and magnetic properties and comparisons 
with experimental data are drawn in section |Vj while 
section [VI| is devoted to conclusions and perspectives. 



II. METHOD 

DFT calculations were performed with the PWSCF 
code of the Quantum ESPRESSO suited in a plane- 
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wave/ultrasoft-pseudopotential approach, adopting the 
Generalised Gradient Approximation (GGA) from Ref!^ 
for the exchange-correlation functional. We used a ki- 
netic energy/charge cut-off of 30/300 Ry. Edges were 
simulated within a super-cell geometry using a vacuum 
layer of 9.5 A between two edges and of 10.6 A between 
two graphene planes. L was fixed according to the lattice- 
constant of graphene (2.46 A) and atomic positions were 
allowed to fully relax. In the case of armchair edges, the 
width of the ribbons was about 22.4 A and for zigzag 
ribbons about 28.9 A. Electronic integrations were done 
using equidistant, consistent grids along the periodic di- 
rection of the ribbon corresponding to 12 k-points in the 
smallest unit cell. Denser grids have been occasionally 
used in the presence of metallic/magnetic states. We lim- 
ited our study to the thermodynamic stability of neutral 
systems, and as such kinetics and charge effects are not 
included. In fact, although graphene edges and GNR can 
be obtained experimentally in charged atmospheres, their 
technological interest resides in their electronic proper- 
ties in the neutral state, i.e. at a Fermi level as close as 
possible to the Dirac point of the graphene band struc- 
ture. We argue that GNRs and graphene edges obtained 
at high doping level might be extremely reactive when 
brought back to neutrality, and recombine to form more 
stable edge structures. 

III. OXYGEN AND CARBON DIOXIDE 

As in RefpSl, we calculate the zero-temperature edge 
formation-energy per unit length as follows: 

£o ^ k ( E " bb ~ NcEblk - *f Eo ) ' (1) 

where E rlbb , E blk and Eo 2 are the total energy of the rib- 
bon super-cell, of one atom in "bulk" graphene and of the 
isolated O2 molecule. Ac (No) are the number of carbon 
(oxygen) atoms in the super-cell, while L the length of 
the unit cell (see Fig. [I]). £o 2 can be used to determine 
the stability of different structures as a function of the 
experimental conditions^. In presence of molecular O2 
gas, at a given chemical potential po 2 , the relative sta- 
bility is obtained by comparing Go 2 = £o 2 ~ /°oMo 2 /2, 
where po = No/(2L). At the absolute temperature T 
and for a partial O2 pressure 

Po 2 = H°(T) - H°(0) - TS°(T) + fc s Tln (j^j , (2) 

where H° (S°) is the enthalpy (entropy) at the pres- 
sure P° — 1 bar obtained from RefH^. In presence of 
monoatomic-oxygen gas one should use Go = £0 ~ Po^Oi 
where £0 — £o 2 ~Po x 2.86 eV, i.e. the binding energy per 
atom of the O2 molecule in DFT-GGA. The most stable 
structures are then obtained, as in Ref.^, by comparing 
their respective Gibbs free energy Go 2 as function of the 
(molecular) oxygen chemical potential. G is linear in /1, 



the slope being determined by po ■ For a given value of p 
the most stable structure is the one with the lowest value 
of G, thus, by increasing p (i.e. going to an environment 
richer in oxygen) the favorable structures will be those 
with higher oxygen-density po- 

We report in Fig.[T]a pictorial representation of the dif- 
ferent GNR edges we have considered in this study, along 
with their characteristic periodicity L. In the left top line 
we show the clean zigzag and armchair GNRg22. The 
right top line (second line from top) represent the zigzag 
(armchair) oxygen-covered nanoribbons included in this 
work. Structures are chosen as to saturate each oxygen 
atom with two single bonds or a double one, while other 
possible chemical bonding configurations are much less 
stable and are not included in the results. The edges are 
ordered with respect to increasing oxygen density and, at 
equal density, with respect to the energetics, fr om m ost 
favorable to least favorable. Non-aromatic edges^SH are 
indicated with a tilde (see section [V]) . 

The density and zero-temperature edge-formation en- 
ergy are reported in Table [I] The free-energy G stabil- 
ity diagram is reported in Fig. [2j We note that three 
zigzag GNRs, named z Q 2, z D 4, and z D 5, and three armchair 
GNRs, named a n, a D i2, and a Q i3, have negative free 
energies, i.e. graphene should spontaneously break and 
form such ribbons if the experimental conditions allow 
to overcome the reaction barriers, as for example when 
nanoparticles, while sliding on graphene layers, are able 
to cut edges along their patrpESl w e note that those ex- 
periments are performed in hydrogen atmosphere; in the 
case of oxygen atmosphere the reaction is more exother- 
mic and could therefore be less controllable. 

The thermodynamic stability diagram (Fig. 2) indi- 
cates that the armchair a i2 and the zigzag z D 4 are the 
most stable configurations at low and high concentrations 
of atomic or molecular oxygen, respectively. However, 
the former one is more stable for nearly all negative val- 
ues of the O2 chemical potential, which suggests that the 
armchair edge should be the only thermodynamically sta- 
ble one in ordinary experimental conditions going from 
ultra- vacuum to atmospheric concentration of molecular 
oxygen. 

On the other hand, it is even more important to con- 
sider the presence of CO2 molecules, either from the at- 
mosphere or from the oxidation of edge carbon atoms. 
We can define, analogously to Eq. [T] an edge formation 
energy per unit length with respect to carbon dioxide, 
which reads: 

(3) 

By simple calculations one can easily show that 

£co 2 = £o 2 + PoAE C o 2 , (4) 

where po is, as above, the oxygen density per unit 
length, and AE C o 2 = (E blk + Eo 2 ~ E C o 2 ) is the for- 
mation energy of molecular CO2 from graphene and oxy- 
gen molecules, and equals -3.76 eV in DFT-GGA. In 
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other words, the chemical reaction C blk + O2 — !• CO2 
is strongly exothermic, as well known from the fact that 
ordinary coal burns in oxygen once ignited. As a con- 
sequence, the molecular-oxygen thermodynamic stability 
scale equally holds for carbon dioxide, with a horizon- 
tal shift of AEco 2 - Conversely, the molecular oxygen 
thermodynamic stability scale must include a vertical 
line at ^o 2 = AEco 2 = —3.76 eV: at higher oxygen 
densities no edge decoration can be thermodynamically 
more stable than CO2. This, again, implies that only 
the armchair a Q i3 edge is realistically observable at ther- 
modynamic equilibrium and oxygen concentration lower 
than the one which would induce spontaneous ignition of 
graphene with oxygen to produce CO2. 



would not spontaneously break to form edges in a wa- 
ter atmosphere. On the other hand, any freshly-cut edge 
would be substantially stabilized in the presence of water 
by the saturation of its dangling bonds with O/H/OH 
groups. The other result is that z hi is the most sta- 
ble configuration at any chemical potential. This means 
that, in principle, the presence of oxygen favors the arm- 
chair edges, while the presence of water favors the zigzag 
ones. In the case of ammonia, free energies are always 
positive, and the z n hi zigzag GNR is the most stable one, 
in correspondence to z a hi in presence of water. 



V. DISCUSSION 



IV. WATER AND AMMONIA 

Besides oxygen, we considered GNR edges saturated 
with water or ammonia molecules. In these cases we re- 
stricted our study to O/H/N contents as fixed by the 
respective molecular stoichiometries, which corresponds 
to the description of experiments carried out in pure wa- 
ter or ammonia atmospheres. In principle, the full study 
of other stoichiometries would allow a more complete de- 
scription of the thermodynamic properties as function of 
two independent chemical potentials, i.e. O and H (N and 
H) in the case of water (ammonia). However, the num- 
ber of different possible edges is formidable, which would 
make these calculations extremely demanding, virtually 
impossible to carry out "by hand" , and likely needing 
random-search optimizations 41 . Moreover, we only con- 
sidered aromatic edge configurations^, with a few ex- 
ceptions in the case of water. The zero-temperature edge 
formation-energies per length reads as follows: 

£h 2 o = ^ (E ribb - N c E blk - N U20 E U2 o) (5) 

and 

£nh 3 = ^ ^E nbb _ N c E blk — A/NH3-ENH3) (6) 

in the case of water and ammonia, respectively. The 
restriction on aromatic edges implies that zigzag edges 
respect the 2-1-1 bond order rule fo r the dangling carbon 
bond that, as shown in RefsP^22l j s the only aromatic 
configuration for zigzag GNRs. In the case of water, sin- 
gle dangling bonds are thus saturated by H atoms or 
OH groups, while double bonds are saturated either by 
oxygen atoms, or by H/OH-containing pairs (see Fig. [I] 
bottom row). In the case of ammonia, we considered all 
water-decorated structures, and replaced each O atom 
(OH group), by a NH (NH2) group. The formation en- 
ergies are summarized in table |TJ the corresponding ther- 
modynamic stability plots are shown in Fig.|3]and Fig. [4] 
In the case of water, two main results are evident: the 
free energy is always positive, meaning that graphene 



We have shown in our previous works on hydrogen- 
decorated edgesSSiand GNRs^that a number of different 
edge structures with different electronic/magnetic prop- 
erties compete for thermodynamic stability at different 
hydrogen concentrations, and that they should be con- 
sidered in order to provide results of practical use to 
experimentalists. In particular, we showed that an im- 
portant hint on the structural and chemical stability of 
GNR edges comes from the well-known organic chem- 
istry structural properties of polycyclic aromatic hydro- 
carbons (PAHs) , large molecules consisting of fused car- 
bon rings. According to the so-called Clar's rul e 1 42 * 43 !, the 
most stable structures for PAHs are the ones maximizing 
the number of Clar sextets. 

In bulk graphene all 7r-electrons take part in the for- 
mation of Clar sextets. In this arrangement, one out of 
3 hexagons is a Clar sextet. We define, as in RefsP^U, 
aromatic edges as those allowing for the same 1/3 den- 
sity of benzenoid rings as bulk graphene, and in which 
every carbon atom is fully saturated (it has 4 covalent 
bonds). Conversely, in non-aromatic edges the 1/3 bulk- 
graphene-density of benzenoid rings is incompatible with 
a full saturation of all carbon dangling bonds. These 
edges are indicated with a ~ in their label. This den- 
sity can only be recovered at the price of having carbon 
atoms with more or less than 4 covalent bonds. In the 
presence of magnetism such over/under-bonding can be 
interpreted as radicals, as in the edge labelled & q (see 
Fig. [5J panel B, and discussion in RefP^). 

As shown in Fig.[T]and in Table[j] we find that, except 
for z Q 2 and analogously to the case of hydrogen-saturated 
edges, the energetically most stable ones are aromatic. 
In the case of water, the z /j4 edge is the only relatively 
stable non-aromatic configuration. In Fig. [5] we report 
the electronic band structure, close to the Fermi level, of 
the most stable and interesting edges. The gray-shaded 
areas indicate the edge-projected bulk bands. The full 
lines correspond to GNR electronic states. Bands outside 
the gray-shaded area are electronic states exponentially 
localized at the edges (edge states'). 

In line with the results of RefP^, in Table []] and in 
Fig. [5] we show that non-aromatic edges invariably dis- 
play edge states at the Fermi level that can split up and 
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thus induce a magnetic character of the edge. In absence 
of such splitting those states are only partially occupied, 
giving rise to metallic behavior. On the contrary aro- 
matic edges either do not possess edge states (as shown 
in Fig. [5| panel A) or such states are fully empty or oc- 
cupied (as acP^) ■ 



VI. CONCLUSIONS 

In conclusion, our ab initio DFT calculations indicate 
that neutral graphene edges, in the presence of molecular 
species typical of experimental conditions, are stabilized 
by decoration with double-bonding O /NH groups or sin- 
gle bonding OH/H/NH2 groups. In contrast to the case 
of hydrogen saturation, for each molecular species we ob- 
serve only very few possible edges in the thermodynam- 
ical stability diagram, namely an armchair edge in the 
case of oxygen, and one zigzag edge in the case of wa- 
ter and in the case of ammonia. This observation might 
help experimentalists in discriminating between zigzag 
and armchair edges and their mixtures according to the 
chemical composition of the atmosphere in which GNRs 
are obtained. 

Calculations were done at IDRIS (project n° 081202). 



TABLE I: Free energy £ (defined as in Eqs. {l}, (5}, (6}), 
oxygen density (po = No/(2L)) and nitrogen density (pN = 
-/Vn/(2L)) for all the studied edges. The z 2 edge, non- 
magnetic within DFT-GGA calculations, is found to be mag- 
netic when hybrid DFT functionals are usecP^. 
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FIG. 1: Scheme and labelling of the oxygen- (right top line and second line from top), water- (third and fourth lines from 
top), and ammonia-terminated (bottom line) edges of a graphene nanoribbon considered in this work. Hydrogen, oxygen and 
nitrogen atoms are the small circles, while carbon atoms are not explicitly represented; dashed lines indicate hydrogen bonds. 
The structures are periodic along the ribbon-edge with periodicity L, and are ordered with respect to increasing O and N 
density and, at equal density, to increasing formation energy. A tilde over an edge label indicates that that structure is not an 
aromatic one. 
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FIG. 2: Free energies versus chemical potential for the five 
most stable oxygen-terminated edges, and of the clean one. 
The alternative bottom (top) axes show the pressure, in bar, 
of molecular O2 (carbon dioxide) gas corresponding to the 
chemical potentials at T = 100, 300, and 600 K. 
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FIG. 3: Free energies versus chemical potential for the eight 
most stable edges after exposure to water, including the clean 
one. The bottom axes show the pressure, in bar, of molecular 
H2O vapor corresponding to the chemical potentials at T = 
100, 300, and 600 K. 
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FIG. 4: Formation energies versus chemical potential for the 
four most stable edges after exposure to ammonia, and of 
the clean one. The bottom axes show the pressure, in bar, of 
molecular NH3 vapor corresponding to the chemical potentials 
at T = 100, 300, and 600 K. 
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FIG. 5: A) Scheme and electronic band-structure of the five 
most stable oxygen-, water- and ammonia-terminated edges 
of a graphene nanoribbon. Carbon-carbon bonds are repre- 
sented with the standard notation, while hydrogen, oxygen 
and nitrogen atoms are the small circles. The structures are 
periodic along the ribbon-edge with periodicity L. The gray 
area corresponds to the electronic-bands allowed in "bulk" 
graphene. The dashed line is the Fermi level. The top-right 
inset displays the standard representation of the aromatic car- 
bon ring. B) Three examples of "metallic/magnetic" edges 
(see the definition in the text). z D 2, non-magnetic within 
DFT-GGA calculations, is found to be magnetic when hy- 
brid DFT functionals are usecP^ z ha and 5 6 both magnetic. 
The fi 6 structure is non-aromatic according to our definition, 
since the two small dots represent a diradical, which is the 
only representation compatible with the carbon chemical va- 
lence. 



